clear all
clc
load V;
gama = 1e-3;
V = sparse(V(:,1:10));
[A, B, E, H, N] = readD();
v=orth(randn(size(V)));
Ar = v' * A * v;
Br = sparse(randi([1 10], 1, 20), randi([1 2], 1, 20), rand(20, 1), 10, 2);
Cr = Br';
Hr = sparse(randi([1 10], 1, 200), randi([1 100], 1, 200), rand(200, 1), 10, 100);
Nr = v' * N * v;
C = B';
[Ar,Br,Cr,Nr,Hr] = GSI(A,B,C,gama*N,gama*H,V,Ar,Br,Cr,Nr,Hr);